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Abstract —This paper investigates a three-node amplify-and- 
forward (AF) multiple-input multiple-output (MIMO) relay net¬ 
work, where an autonomous relay harvests power from the source 
information flow and is further helped by an energy flow in 
the form of a wireless power transfer (WPT) at the destination. 
An energy-flow-assisted two-phase relaying scheme is proposed, 
where a source and relay joint optimization is formnlated to 
maximize the rate. By diagonalizing the channel, the problem 
is simplified to a power optimization, where a relay channel 
pairing problem is solved by an ordering operation. The proposed 
algorithm, which iteratively optimizes the relay and source power, 
is shown to converge. Closed-form solutions can be obtained for 
the separate relay and source optimizations. Besides, a two-phase 
relaying without energy flow is also studied. Simulation results 
show that the energy-flow-assisted scheme is beneficial to the 
rate enhancement, if the transmit power of the energy flow is 
adequately larger than that of the information flow. Otherwise, 
the scheme without energy flow wonld be preferable. 

Index Terms —Energy harvesting, multiple-input multiple- 
output (MIMO), relay network, amplify-and-forward (AF). 

I. Introduction 

As a promising technology for energy-constrained wire¬ 
less networks, joint wireless information and power transfer 
(JWIPT) now attracts much attention in the context of relay 
networks. 

Current research on JWIPT in relay networks mainly studies 
single-antenna systems with energy-constrained relays Q-Gl 
or energy harvesting issues in multiple-antenna systems with 
power-supplied relays a, 0. Relaying protocols for wireless- 
powered relays were firstly proposed in m, which focuses on 
3-node single-antenna relay systems by applying the unified 
power splitting (PS) and time switching (TS) frameworks 0. 
More general scenarios of multiple source-destination pairs 
were studied in a. Besides the above works on single-antenna 
systems, a studies a 3-node one-way MIMO relay system 
where a separated energy harvester extracts wireless power 
from signals transmitted by a source and a relay. In a, a 
multiple-antenna relay with power supply is considered. Al¬ 
though a studies a wireless-powered multiple-antenna relay, 
all other nodes are assumed single-antenna, and the processing 
matrix at the relay is not optimized. Different from the 
above research, this paper investigates a one-way amplify-and- 
forward (AE) multiple-input multiple-output (MIMO) relay 
network with a wireless-powered relay. 

As shown in Fig. [T] we consider the scenario where there is 
no direct link between S and D due to barriers (which causes 
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Fig. 1. JWIET relay network. The destination, relay, and source are 
designated as D, R, and S, respectively. 

huge shadow fading), such that the transmission between those 
two nodes has to rely on a wireless-powered relay R. Different 
from the previous works, we also consider the simultaneous 
transmission of power (WPT) and information from D and S, 
respectively. In order to efficiently utilize the energy flow (i.e. 
WPT) but not increase the timeslot consumption, a two-phase 
energy-flow-assisted relaying scheme is proposed. To make the 
formulated rate maximization optimization tractable, it is sim¬ 
plified to a power optimization by performing a channel diag- 
onalization based on a harvested-power-maximization power- 
leakage-minimization (HPM-PLM) strategy. Power allocation 
at R and S are optimized based on an alternating optimization 
(AO). Channel pairing issues introduced in the relay power 
optimization are solved by an ordering operation. Closed-form 
solutions can be achieved in the separate relay and source 
power optimizations. While the energy flow provides the relay 
with an additional source of energy to amplify and forward the 
information flow, the latter is now subject to the interference 
from the energy flow. Hence part of the energy harvested at 
the relay is consumed to amplify and forward the interference, 
which reduces the power usage effectiveness. An alternative 
strategy would be to simply rely on a two-phase relaying 
without (the support of) energy flow. Simulation results in¬ 
dicate that the energy flow-assisted strategy is beneficial to 
the rate enhancement if the transmit power of the destination 
is adequately larger than that of the source. Otherwise, the 
two-phase relaying without energy flow would be preferable. 

The remainder of this paper is organized as follows. The 
system model of the energy-flow-assisted two-phase relaying 
is formulated in Section [II] Section |In] then performs channel 
diagonalization. Section|IV]discusses the joint relay and source 
power optimization. Section lY] elaborates on the two-phase 
relaying without energy flow scheme. Section [Vl] evaluates the 
performance of the schemes. Finally, conclusions are drawn in 
Section IVIII 

Notations: In this paper, matrices and vectors are in bold 
capital and bold lower cases, respectively. The notations (A)^, 







(A)*, (A)'^, Tr{A}, det (A), Ai (A) and [A]i represent the 
transpose, conjugate, conjugate transpose, trace, determinant, 
the ith eigenvalue and the ith column of a matrix A, re¬ 
spectively. The notation A ^ 0 means that A is positive- 
semidefinite, and 7r(a) and ||a|| denote the permutation and 
2-norm of a, respectively. When ^ and are used, top cases 
or bottom cases in the two notations hold simultaneously. 

II. System Model and Problem Formulation 

In Fig. [T] each node is equipped with r antennas. The D- 
R, S-R, and R-D channels are respectively designated as 
Hr.d e e and 11d,r S which 

are independent and identically distributed (i.i.d.) Rayleigh 
flat fading channels, and all the channel matrices are full- 
rank. Due to channel reciprocity, Global 

CSIT is available at each node. The relay exploits PS scheme 
ID for simultaneous energy harvesting (EH) and information 
detecting (ID). At each antenna of the relay, a fraction of the 
received power, denoted as the PS ratio pm for m = 1,..., r, 
is conveyed to the EH receiver. In this paper, uniform PS is 
assumed, i.e. pi,..., = p. The noise at the ID receiver 

(at R) and D are respectively denoted by n/{ ~ CAA(0, cr^I) 
and no ~ CAA(0, cr^I), while the effect of noise at the EH 
receiver is small and neglected ii-ii. 

In phase 1, the received signal at the EH receiver is given 
by yR,EH = {Ur^d^d + Hfi,sxs), where and xs 
are precoded signals from D and S. Assuming an RE-to- 
DC conversion efficiency of 1, the harvested power equals 
TrlpH/i^DQ^Hf ^-f pHfl_sQsHf g}, where = 

£{xo^o} Qs = respectively. Meanwhile, 

the baseband signal input to the ID receiver for forwarding 
is given by = (1 - p)^/^ (Ur^d^d + 'H.r,s^s) + ^r. 

In phase 2, the information received at D is given by 

yD=(l —(Hi{_sxs + Hi{_£)X£))-f n^-l-n£),(l) 

where = Hupirni} and F denotes the relay process¬ 
ing matrix. With perfect channel state information, the self¬ 
interference in yo, i.e. the term related to xd, can be canceled, 
but some power at the relay is consumed to forward this self¬ 
interference. To maximize the achievable rate, an optimization 
problem can be formulated as 

PI: max i logdet (l-f (l-p)H, 5 .itrHfl,sQsHf gF^- 
Qd.Qs.F ^ 

(2a) 

s.t. Tr {(1 -p) (FHfl.sQsHf_sF^ + FHfl,z?QDHf^^F^) 
+alFF^}=pTT + ,(2b) 

Ti'IQd} < Po , Qd h 0 ,Tr{Qs} < Ps ,Qs hO, (2c) 

where the coefficient of 1/2 in (|2a|l results from the half¬ 
duplex transmission, W = a‘^Flo,RFF^Flj^ + cr^I, and 
(l2bl) implies that all the harvested power at the relay is used 
for forwarding. For simplicity, p is not optimized, but an 
exhaustive search is conducted to find the best p. Although 
the simultaneous transmission in phase 1 is similar to the two- 
way MIMO relaying |[8|, optimizing matrices directly as in IS] 


is intractable due to the mutual information criterion ( l2a] i and 
the constraint (l2bl) . To solve the problem, an iterative algorithm 
based on channel diagonalization is then proposed. 

HI. Channel Diagonalization 

To simplify the design problem, this section decomposes the 
forwarding channel in phase 2 and jointly decomposes 

the S-R effective channel = H/j and the D-R 

channel H/jpj in phase 1 based on the HPM-PLM strategy, 
such that (l2al) and (l2bl) can be diagonalized and problem PI 
reduces to a power optimization problem. 

A. Structure of Relay Matrix 

As an unique forwarding channel, FLo,r is decomposed 
as its singular value decomposition (SVD) XJo,r'^d,rYo r- 
Then, with the SVD of FIr^s = ^r,S^r,S^r^S’ ^PPlyi^g 
the matrix inversion lemma to ( l2al ) yields 

G = i log det f I + (I - (1+ 

2 V CTn 

, (3) 

where the matrix between the two equals a positive 

semidefinite matrix Uf ^(I-t-HD,_RFF^Hg ■ 

HupjFU/j g. Hence, the matrix in det(-) is positive-definite. 
According to Hadamard’s inequality IH, (O is maximized pro¬ 
vided Uf 5 F^Vr,,RS|,_^Vg^FU/i,s is diagonal. Hence, 
F = 'Vo,r^f^r s for C which means that the 

relay couples a given receive eigenmode of U_r,s with a given 
transmit eigenmode of V^pj with an amplification factor 
given by the corresponding diagonal entry of With the 
decomposed F, (l2al l can be diagonalized. 

B. Maximize Harvested Power and Minimize Power Leakage 

The HPM-PLM strategy is proposed to enhance the power 
usage effectiveness at R, since the energy flow harvested at 
the EH receiver is used to forward not only the information 
but also the energy flow leaking into the ID receiver, as shown 
in (Hbl. Performing eigenvalue decompositions (EVD), Qd = 
VDSgVg and HD,DQDHf D=UF,DS?j,DUf,D- Given 
the SVD of Hd,d = Vg rearranging ( l2b] i yields 

Tr |(1 — g + cr^SgSD — (4) 

= Tr {(pi - (1 - p)Sf Ed) (ug,gVg^«SD.RUg_^VD- 

(5) 

= Tr{(pI-(l-p)SfSD)(uI,sUD.DSg,^Ug^Ufl,s)}.(6) 

Eq. (|6]l highlights the difference between the power of the 
energy flow harvested at the EH receiver and the energy 
flow leaking into the ID receiver. Thus, a strategy can be 
proposed to maximize the harvested energy and minimize the 
power leakage. To maximize the transfer of energy from D 
to R, rank-one transmission should be exploited at D, i.e. 
Qd = PdIFJo rU^xIFXo r]!^^^, where [U^p^]max is the right 



singular vector (RSV) corresponding to the maximum singular 
value of H_r,_d, such that pTr{H_R,_DQ£)Hf = 

pPD^DM,m^x (see Proposition 1 in CqI for proof). To min¬ 
imize the power leakage, the power leaking into the ID 
receiver should be paired with the minimum amplification 
coefficient, such that the power of the retransmitted leakage 
equals (1 - p)Xf,mmPD^D,R,max, where A/,min denotes the 
minimum diagonal entry of With the above strategy, 

(l6ll should equal (p - (1 - p)Xf^tnia)PD>>D,R,nmx, which is 
shown to be an upper bound of © by applying Lemma 
ii.i in mi]. To make © equal to the upper bound, in ©, 
= P.. This specific P^r permutates the unique 
non-zero diagonal entry PnXD.R.max in the diagonal matrix 
Sd to the same position as 

p — (1 — p)A/,niin in pi—(1 — p)S|^Si?. By this means, © 
achieves the upper bound, and (l2Ft is diagonalized. 

In summary, (|2a]l is diagonalized with the decomposed 
a-D,R and the structure of F, and (l2Ft is diagonalized 
with the HPM-PLM strategy, i.e. the rank-one Qd and 
Ufl,s = VJ, pPl. Since Qs = where = 

Tr{Qs} = l|he.™|pArj,s.™ < Ps, 

where he,m = [He]m and \R^s,m denotes the mth diago¬ 
nal entry of g. Hence, problem PI reduces to a power 
optimization. Since F = V the RSV of 
F (i.e. g) matches the left singular vectors (LSV) 

of and a permutation of the LSV of /{. Accord¬ 

ing t( 

P.S 

where 

because of the same (but permutated) LSV of H/{_s and 
u, the channel power gains of the effective S-R and D- 
R channels in phase 1 are overlapped at the ID receiver, i.e. 
(l-p)(Sls + P^^i),riUS,flQ£'Ui5.flPj)’ where the diag¬ 
onal entries of g and (1 — p)P 7 rS |5 
are respectively denoted as Xr^s — [Afl,s,i, ■. ■, 
and j3 = [/3i,...,with the unique non-zero /?„ = 
(1 — p)PDA£)_r{_niax = c. Although the retransmitted energy 
flow can be canceled at D (i.e. P^S|j rQd"U*d 
in £{yDyD} is deleted), the overlapped channel power gains 
in phase 1 still impact the rate, because they are amplified 
and transmitted. As shown in the diagonalized relay power 
constraint (l2b] l (i.e. the following dTell). the diagonal entries of 
denoted as A/ = [A/p,..., A/,r]^, are non-uniformly 
weighted by the entries of the overlapped channel power 
gains i.e. (1 - p)XR^s,m + /3m for m = l,...,r. Thus, 
the pairings of the diagonal entries of ^ (related to the 
forwarding eigenmodes in phase 2) and the overlapped channel 
power gains in phase 1 affect the optimization of A/ and 
thereby the rate. Additionally, in phase 1, the value of each 
(1 — p)Xpi^s,m + /3m is affected by (which determines the 
pairing of each Ai{,s.m and /3m)- 

IV. Joint Power Allocation Optimization 

To make further calculation and analysis tractable, we then 
focus on the achievable rate at high receive SNR. Substituting 


I O, £{yDy^) = VD,R[T.l^pT.pT.f{i-p){T.l^g + 
1,r^I,r^d\5*d,rPI) + + l)a2]Ug 

“ Ujj/{Q liUij is diagonal. In the above f{y_Dy^}, 


the previous channel decompositions into problem PI, the 
original problem can be reformulated as 


m ■ 1 ( (1 “ P)A_R,S,mA/,mAD.R,i 

P2: mm - > log -——-T-V 

A/,Ab,s (1 -I- A/,mAD,fl,mj 

s.t. A/p, A/p,..., A/_r > 0 , 

Aij.sp, A_r,5P, . . . , Xr^Sp > 0 , 

|^A_R,5,m < Ps, 


Tr{Qs} = E 


m—1 


(7a) 

(7b) 

(7c) 

(7d) 


r 

((1 - P)XR,S,m 

m—1 

pPoXo ,R,max ; 


r 

+ pXR^S,m + 

m—1 

(7e) 


where /3m is constrained by /3m = c, if m = index(A/^min) 
(where index(A/^tnin) returns the index of A/^min); otherwise, 
/3m = 0. Problem P2 is not convex due to the non-affine (iTel) . 
Then, problem P2 is solved using an AO. 


A. Relay Optimization with Fixed Source Power Allocation 

With given XR^s,m, the power optimization problem at the 
relay is formulated as 


P3(a): max log [ 

m=l \ 

S.t. (fTbl i and (iTel) . 


(1 - P)XR,S,mXf ,m Ad ,_R 

O’™ (1 + Xf^mXn 



( 8 ) 


The challenge in solving P3(a) is that c of Pm is required 
to be paired with A/^min, but the position of A/^min in A/ is 
unknown before solving the problem; the rate is affected by 
the pairings of the elements of Xd,r = [Arj^ijp,..., Xp^RpY 
(i.e. the diagonal entries of ^ related to the forwarding 
eigenmodes in phase 2) and Xr^s (related to in phase 

1), even if the constraint on c is relaxed and /3m is fixed. To 
avoid the high complexity of searching the best pairings for 
P3(a), we then reveal that the pairing issues can be solved by 
ordering operations. Firstly, we relax the constraint on c, i.e. 
P^n- becomes an arbitrary permutation matrix P^r, and assume 
columns of V d^r are arranged in certain orders, such that the 
elements of Xp R and Xr^s paired in certain ways and c 
is paired with a certain XR^s,m- Problem P3(a) then becomes 
convex regardless of the pairing issues. By analyzing Karush- 
Kuhn-Tucker (KKT) conditions, a closed-form solution can be 
obtained by 


•^/,m — 

^XD,R,r. 


^ ^D,R,m V*XD,R,m — p)Xr^S, m +Cr'^+ 


,(9) 


where v* denotes the Lagrange multiplier for the constraint 
Gi and is greater than 0. It can be calculated by solving 

Em=l '^/,m P)Xr^B ,mP '^n P /3m^ = pPoXp ^R^max P 

Em=i P^R,S,m with bisection. By using ©, two lemmas 
are revealed. For notational simplicity, it is defined that 









^ Im = 


Z'm — (1 P)^R,S^m “t“ ' 

(1 - p)\R,S,m + Cr^, 1 = [/i, . . . , 

Lemma 1: Suppose that the elements in 7ri(z) are arranged 
in the same order as another permutation 712 ( 2 ) except that Zi 
and Zj (where Zi < Zj for i < j) in 7ri(z) are swapped in 
772(2), i.e., Zi and Zj in 7 ri(z) are paired with Xd,r,p and 
\D,R,q (where ^D.,R,p < ^D,R,q for p < q), respectively, 
while Zj and Zi in 772(2) are paired with Xd,r,p and XD.R,q, 
respectively. Then, the value of the objective function ® with 
771 (z) is no less than that with 772 ( 2 ). 

Proof: To prove this lemma, 7ri(z) and 772 ( 2 ) are respec¬ 
tively substituted into (|9l) with Xd,r to calculate and compare 
the values of I®. The lemma is finally justified by scaling 
Inequalities. See our extended version ifTSIl for details. ■ 
Lemma [T] addresses the pairings of the transmit eigenmodes 
of F (i.e. 'VDjj) and the overlapped channel power gains. 
With hxed pairings of Xii,s,m and j3m for m = 1... r, the 
values of the entries of the overlapped channel power gains, 
i.e. (1 — p)Xii^s,m + Pm for m = 1.. .r, are fixed. Lemma 
[T] reveals that, for two pairs of the transmit eigenmodes of 
V o.R and the overlapped channel power gains (while other 
pairings are fixed), the strongest eigenmode of V o,r and 
the strongest overlapped channel power gain should be paired 
together. The following Lemma |2] addresses the pairings of 
the channel power gains of the effective S-R channel and the 
non-zero channel power gain of the effective D-R channels. 


respectively. Hence, the permutation matrix = I = P^. 

B. Source Optimization with Fixed Relay Power Allocation 
According to Proposition [T] Xf,m are arranged in a decreas¬ 
ing order, and index(A/^tnin) = r. Thus, the source power 
optimization problem is formulated as 

1 / P')Xr S,mXf mXo R m V,,., , 

P3(b): mm - > log ^ ^ ' ',’ - 7 ^ KlOa) 

Ail,s \ CTn (1 + Xf^mXD,R,m) J 

S.t. 0 < Xr^sp < Xr^s,2 < • ■ • < XR^s,r , (10b) 

(iTdl i and (17^ . 

Problem P3(b) is convex and can be solved by an optimization 
solver. Analytical solutions are still attractive due to its low 
complexity. The challenge in deriving a closed-form solution 
is the ordering constraint in (llObl l. We hnd that when the 
ordering constraint in (llObl l is relaxed, the output X^ g ^ can 
still be in an increasing order if XR^s.m in (iTdl i are uniformly 
weighted (otherwise, the Lagrange multiplier for (ITdl i would 
be non-uniformly weighted in the derived closed-form solu¬ 
tion, which may violate the ordering of Xr^s)- Therefore, 
we reformulate problem P3(b) as P3(c) by only replacing 
constraints (llObl l and (ITdl i with Xr^s,i, Xr^s, 2 : ■ • ■, XR^s,r > 0 
and EL=imax{|| he,?7i |P}A_R,s.m < Ps, respectively. The 
max{||he,m|P} is denoted as hf in subsequent parts. Prob¬ 
lem P3(c) is convex, and KKT conditions are given by 


i.e. the pairings of XR^s,m and c for m = 1... r. It is shown 
that, for two channel power gains in XR^s,m, the strongest 

^*R,S,m > 0 I 
r 

m = 1,..., r 

(11) 

XR,s,m should be paired with c. 

Lemma 2: Assume two permutations 7 ri(l) and 772(1). In 

^ ^ ^R,S,m 
m—1 


(12) 


771 (1), positions of p and Ij follows that min{Zi + c,lj} and 
max{(i-|-c, Ij} are respectively paired with Xd^rp and Xr^rj 
(where i < j, p < Ij and Xd^rp < XD,R,j)- In 772 ( 1 ), positions 
of li and Ij follows that inm{li,lj -f c} and ma,x{li,lj + c} 
are respectively paired with XD,R,i and Xd^rj- Other pairings 
between XD.R,m and Im (for m ^ i,j) in 7ri(l) are the same 
as 772 ( 1 ). Then, the objective function ® with c paired with 
Ij yields a higher value than that with c paired with p. 

Proof: Based on the conclusion of Lemma [T] Lemma |2] 
is proved. The idea to prove this lemma is similar to Lemma 
m See na for details. ■ 

Proposition 1: When the elements in Xr^s and Xd,r are 
arranged in increasing orders and the non-zero pm is paired 
with the maximum XR^s,m, the value of the objective function 
® in Problem P3(a) is maximized and the optimal AJ ^ are 
arranged in a decreasing order. 

Proof: Applying Lemma [T] and Lemma ^ the orderings 
of Xr^s and Xjj r and the pairing of Pm and XR^g^m in the 
proposition are proved through induction. Based on the above 
orderings and pairing, the ordering of Aj ^ can be easily 
proved. See the Proposition [T] in 1(121 for details. ■ 

Proposition [T] illustrates that the constraint on Pm (i-C- c 
is paired with A/^min) can be safely relaxed. Following the 
ordering operation in Proposition [T] entries (p — (1 — p)A/^tnin) 
and PdXd ijmax are at lower-right corners of matrices pi — 
(1 -p)SfS;. and in ®, 


^ ^ (•^/,m(l p) P)^R,S,m ~ ^ ^ (fn + Pm)Xf^, 


m—1 


m—1 

+pPdXd 

,R,ma\ ) 

llm > 0 , 

72* >0, 

^2 f X! ^R,S,m - . 2 ^ 0 ’ 

\m=l "■e.max / 


(13) 

(14) 

(15) 

(16) 

(17) 


and 


A 


- - + 72 + P* (A/.m(l - P) - P) = 0, (18) 

R,S,m 

where 7 ^ 72 , and p,* denote the optimal Lagrange multipli¬ 

ers. Eq. (ITO . (fl® . and dTSll reveal that 7 ^ ^ = 0. If 72* = 0 , 
according to (ITSl l. it is obtained that 

= p*{Xpm{l-p)-p) ’ 

where p* is obtained by solving r/p* = pPDXD,R,max — 
+ Pm)Xf,m- Since yRgm > 0 Vto and p* also 
conforms to (ITZl i. ([T9l l is obtained provided 


A/.m(l - p) - P ^ 0 , Vm 

1 < Ps/h e,max 


( 20 ) 












Algorithm 1 Two-phase relaying with P3(b) 


Initialize A 
repeat 


(0) 


and g 


Update 


by calculating (|9l); 


Update A^g^^ by solving P3(b); 

K ■<— K + 1; 


until 




(k) r(K) 


< e 



Fig. 2. Average rate as a function of PS ratio with d^Rldos = 0.65. 


Algorithm 2 Two-phase relaying with P3(c) 


Initialize and A^g 


repeat 

Update by calculating (|9l); 

if ( l20b is satisfied then 

Update A^g^^ by calculating ( fT^ : 
else 

Update A^ 5 ^^ by calculating (|2TI) : 

K <— K + 1; 


until 




(k) T(k) 


< e 


the previous relaying scheme, due to the absence of the 
energy flow Qu, the S-R channel can be decomposed by 
SVD. Recall that the SVD of g, 

where = diag{X[i^s,i, ■ ■ ■, >^R,S,r}’ *e EVD of 

H«,s 9 'sH |5 = where = 

diag{X'j^s,v-,^'R,S,r}- Therefore, Qj, = V/j^sS'gVf 5 , 
where = diag{X'g-^^,...,X'g^^}, U'^ c, = IJr^s, and 
'^'r,s = '^s'^R,s- For simplicity, we only focus on the case 
of uniform source power allocation, i.e. X'g ^ = Ps/r Mm. 
Assuming high receive SNR and omit the coefficient 1/2, the 
power optimization of problem P4 is reformulated as 


is satisfied. On the other hand, if A 3 > 0, the optimal X^ g ^ 
is achieved by 

where 72 and /i* can be obtained by solving the non-linear 
system composed of (fl^ and (foT l. 

As a summary, the proposed AO-based joint optimization 
algorithms are outlined in Algorithms [U and 121 where the 
objective function (Gill is denoted as C{Xf, Xr^s)- Since the 
separated optimization problems at the relay and the source 
(i.e. P3(b) or P3(c)) are convex problems with strictly convex 
objective functions, value of C{Xf, Xr^s) monotonically de¬ 
creases with each iteration. Besides, the objective function (l7al) 
is lower-bounded. Thus, the two algorithms finally converge. 

V. Two-Phase Relaying without Energy Flow 


T,, . , I (1 P)^R,S,m^f,m^D,R 

P5;min—> log - j - - - 7 

"'/ \ <jl[l + X'f^^XD,R,m) 

s.t. Xf i, Aj 2 i ■ ■ • I ^f,r — 0 j (23b) 

r r 

+ pA'^5„|23c) 

m—1 m—1 

where A/ = [A/, A/ The pairings of X'g g ^ and 
XD,R,m for TO = l,...,r can be solved by Lemma [T] 
with Pr = 0 and /3 = 0. Hence, if X'g^ g .^ = Ps/r ■ 
[Ar?,s.min, ■ • ■, Aii,s,max]'^ and Xo.R.m are arranged in an in¬ 
creasing order, the pairing problem is solved. Thus, the optimal 
X'j ^ is obtained by 

[7J* = -xr^ 

2XD.R,m 



Considering that the relay only harvests power from the 
information flow (i.e., only the S-R and the R-D links exist in 
Fig. [B, a two-phase relaying without energy flow is proposed. 
The design problem is formulated as 

P4: max ^logdet fl-|- 
Qg.F' 2 V 

[W']”') (22a) 

s.t. Tr {(1 - p)F'UR,sQ'glilg [F']" + a^F' [F']"} 

= pTr {UR^gQ'gUlg} , (22b) 

Tr{Q's} <Ps,Qs ^0, (22c) 

where W' = a/,liD,RF' [F'f + a/,1, and F' de¬ 

notes the relay processing matrix. Similar to the energy- 
flow-assisted two-phase relaying, channel diagonalization is 
also used to simplify the design problem. Different from 


-t-- 


where [v']* is constrained by (123cl i. 


0 


,(24) 


VI. Simulation Results 

In the simulations, channel matrix FLij is generated by 
FLij = where j represents the small-scale fading. 

The large-scale fading is given by A“| = where dij is 

the distance between nodes i and j and dos — doR + dRs- 
In the simulations, r = 4, c?ds = 10 m, and cr^ = 1 p,W. 

Fig. 12] shows that for a certain doR/dos value, the average 
achievable rate of the proposed two-phase relaying schemes 
firstly increases to a stationary point at PS ratios of 0.72, 0.74, 
and 0.88, respectively. Then, the rate decreases as the PS ratio 
increases. This is because when the PS ratio is small, less 
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Fig. 3. Rate performance under different don/dus ratios. The PS ratio is exhaustively searched among 0.02:0.02:0.98 to maximize the average rate. 


power is harvested for forwarding, which limits the receive 
SNR at D; when the PS ratio is large, less signal power 
remains for the ID receiver at R and the SNR at R decreases. 
Both the above two cases degrade the achievable rate. It is 
also observed that the two-phase relaying with P3(c) is always 
inferior to P3(b). The reason lies in that the modified source 
power constraint niax{|he,mp}Ai{,S.m < Ps in P3(c) 

would limit the received information signal power at the relay, 
and S-R link performance is sacrificed. 

Fig. |3(a)| shows the average rate as a function of don/dos 
ratio, where an exhaustive search is performed to find the 
best PS ratio for each don/dos ratio. It is shown that the 
rate of the scheme without energy flow decreases as the 
relay approaches to D, because R only extracts forwarding 
power from the information flow. Nevertheless, harvesting the 
energy from D, the rate of the energy-flow-assisted scheme 
can increase as doR/dos decreases. However, it is observed 
that when R is close to S (e.g. doR/dos = 0.9 where the 
energy flow slightly contributes to the forwarding power), the 
scheme without energy flow can outperform the energy-flow- 
assisted scheme. This is because to enhance the power usage 
effectiveness, the LSV of is forced to be V^ ^ with 

the HPM-PLM strategy, which makes the source beamforming 
matrix not always unitary. Fig. |3(b)| demonstrates that if R is 
close to D, the S-R link becomes the critical link; thus, a 
lower PS ratio is needed. Otherwise, the R-D link becomes the 
critical link, and a higher PS ratio is needed. Compared with 
the energy-flow-assisted scheme, the scheme without energy 
flow needs higher PS ratios, because its relay power only 
comes from the information flow. 

Fig. |3(c)| studies the asymmetric scenario where the power 
budget at D is increased while that at S is decreased. It 
is observed that the rate of the energy-flow-assisted scheme 
increases as R moves towards S. This is because with ade¬ 
quately large power budget at D, the receive SNR at D can 
still increase, although the harvested power of the energy flow 
at R decreases. Compared with Fig. |3(a)| the energy-flow- 
assisted scheme outperforms (rate-wise) the scheme without 
energy flow at most doR/dos ratios due to the efficient 
utilization of the harvested power. Recall that the energy-flow- 
assisted scheme prefers a high-quality S-R channel, e.g. a 
high transmit SNR at S. Although the transmit SNR at S 
decreases in the case of Fig. |3(c)| the relative difference in 


rate at doR/dDS = 0-^ between the two schemes in Fig. |3(c)| 
is smaller than that in Fig. |3(a)| This illustrates that the rate 
can benefit from the energy-flow-assisted scheme, when the 
power budget at D is adequately larger than that at S. 

VII. Conclusion 

In this paper, we have investigated JWIPT relaying schemes 
in an AF MIMO one-way relay network, where a wireless- 
powered autonomous relay is deployed. Considering possible 
simultaneous transmission of energy and information, we have 
proposed an energy-flow-assisted two-phase relaying and a 
two-phase relaying without energy flow. Simulation results 
reveal that the rate can benefit from the energy-flow-assisted 
scheme if the transmitted power of the energy flow is ade¬ 
quately larger than that of the information flow. Otherwise, 
the scheme without energy flow would be better. 
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